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ABSTRACT 


The  sensitivity  of  a  passive  horizontal-tracking  algorithm  to  variations  in 
input  measurements  is  investigated.  The  algorithm  determines  estimates  for 
depth,  range,  bearing,  horizontal  speed,  course,  and  frequency  for  a  cw 
acoustic  source  moving  with  constant  velocity  at  fixed  depth.  The  receiver  is 
a  horizontal  linear  array  towed  at  a  constant  depth.  Both  source  and  receiver 
move  in  the  upper  portion  of  a  deep  ocean  and  are  separated  by  a  relatively 
short  range.  Dominant  acoustic  signals  are  presumed  to  arrive  along  two 
upper-ocean  ray  paths.  The  algorithm  uses  a  new  combination  of  input  quan¬ 
tities,  including  multipath  information,  Doppler  frequency  shifts,  and  array 
directional  measurements.  Procedures  are  developed  for  analyzing  effects  of 
input-measurement  errors  on  source  localization.  The  robustness  of  the  algor¬ 
ithm  to  small  variations  in  acoustic  measurements  and  environmental  parameters 
is  demonstrated  for  a  variety  of  source-receiver  configurations.  Variance 
estimates  of  position  and  motion  are  obtained  in  terms  of  input-measurement 
variances.  Bounds  on  tracker  performance  are  developed  for  measurements  that 
are  affected  by  noise.  Results  from  the  several  types  of  analyses  corroborate 
the  sensitivity  characteristics  of  the  algorithm. 
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introduction 

Many  procedures  have  been  developed  to  determine  passively  the  motion  of 
an  underwater  acoustic  source.  They  involve  a  variety  of  acoustic  quantities 
measured  at  one  or  more  receivers,  which  typically  consist  of  arrays  of  sen¬ 
sors.  In  familiar  bearings-only  tracking,  direct  measurements  are  made,  at 
several  observation  times,  of  received  signal-arrival  angles.  These  are  sub¬ 
sequently  processed  to  obtain  estimates  of  source  position  and  velocity.1""^ 

A  quantity  which  has  received  considerable  attention  recently  is  the  time 
delay  (or  travel-time  difference)  between  ray  arrivals  at  sensor  pairs. 

For  example,  Ref.  8  describes  how  time  delays  observed  at  a  linear  array  may 
be  converted  into  source-distance  and  direction  information.  Another  measured 
quantity  (or  descriptor)  is  the  received  frequency,  which  is  employed  in 
Doppler  tracking. ^  some  tracking  methods  have  exploited  a  combination  of 
descriptors,  such  as  multipath  time  delays  in  conjunction  with  bearing 
measurements . ^ 

In  a  previous  investigation, 11  effects  of  a  moving  cw  source  on  acoustic 
total  field  received  at  short  ranges,  was  studied.  The  receiver  was  fixed  on 
the  surface,  and  the  total  field  consisted  of  two  bottom-bounce  rays  emanating 
from  the  constant-velocity  source  moving  above  the  SOFAR  axis.  A  key  feature 
was  the  occurence  of  brief,  regularly  spaced  changes  in  total-field  phase 
rate,  referred  to  as  phase  roll.  We  present  here  a  passive  tracking  algorithm 
using  as  one  of  its  descriptors  the  quasi-period  of  phase  rolls.  The 
algorithm  determines  estimates  for  depth,  range,  bearing,  horizontal  speed, 
course,  and  frequency  for  a  cw  source  moving  with  constant  velocity  at  fixed 
depth.  The  receiver  is  a  horizontal  linear  array  which  is  towed  at  fixed 
depth.  Both  ‘source  and  receiver  move  in  the  upper  portion  of  a  deep  ocean  and 
are  separated  by  relatively  short  ranges.  The  horizontal-tracking  algorithm 
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involves  a  new  combination  of  input  information  from  direct  and  surface- 
reflected  ray  paths.  Specifically,  the  descriptors  include  the  multipath 
phase-roll  quasi-period  as  well  as  the  signal-arrival  angle  and  received  fre¬ 
quency  from  the  direct-ray  arrival.  Tn  order  to  show  the  essential  features 
of  our  procedure  in  as  simple  a  context  as  possible,  the  sound  speed  is  taken 
as  uniform  in  the  ocean  region  from  the  surface  to  the  deeper  of  the  source 
and  receiver,  and  the  ocean  surface  is  assumed  horizontal. 

Many  factors  can  contribute  to  measurement  errors  in  underwater  source 
localization  and  tracking.  One  is  the  omnipresence  of  noise,  which  corrupts 
all  measurements  to  some  degree,  and  arises  from  a  variety  of  environmental 
and  systems  causes.  Other  factors  involve  the  receiving  array  of  sensors, 
which  is  subject  to  physical  influences  which  limit  its  performance  in  esti¬ 
mating  descriptors  such  as  signal-arrival  angle.  Examples  are  the  individual 
independent  motions  of  array  sensors12  or  sub-arrays,13  array  element  failure 
and  errors  in  amplitude  and  phase  between  hydrophone  channels,14  and  bending  of 
the  linear  array.16  Furthermore,  aspects  of  systems  software  and  signal  process¬ 
ing  also  contribute  to  descriptor  errors.  For  example,  measurement  and 
processing  times  must  be  small  enough  so  that  source  or  receiver  motion  does  not 
appreciably  affect  descriptor  values.  A  discussion  of  this  characteristic  in 
relation  to  received  frequency  measurements  appears  in  Ref.  9.  Finally,  the 
neglect  of  variations  in  various  environmental  features  involved  in  tracking 
procedures  can  lead  to  errors.  Examples  of  environmental  influences  which 
have  been  investigated  are  sound-speed  or  bottom-slope  variations,16'17  and 
turbulence  in  the  propagation  medium.18  Often,  the  aforementioned  descriptor 
errors  are  either  inherently  random  or  must  be  treated  as  random.  However, 
some  error  influences  can  be  treated  in  a  deterministic  manner.  As  just  one 
example,  Ref.  19  contains  an  analysis  of  non-random  eddy-induced  variations  on 
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source  bearing-angle  determination  by  a  horizontal  linear  array.  It  is 
necessary  as  well  to  analyze  the  effects  on  tracking  procedures  of  these  types 
of  errors,  which  are  not  necessarily  random  but,  in  fact,  may  arise  from 
modeling  simplifications . 

The  principal  purpose  of  this  paper  is  investigation  of  the  sensitivity 
of  the  aforementioned  horizontal-tracking  algorithm  to  variations  in  inputs, 
including  both  acoustic  descriptors  and  other  algorithm  parameters.  The 
analysis  of  effects  of  input  measurement  errors  on  output  predictions  for 
source-localization  quantities  can  typically  be  divided  into  two  stages.  The 
first  is  determination  of  input  variations  resulting  from  the  specific  error 
source,  and  the  second  is  the  study  of  input  variations  on  various  outputs 
using  the  algorithm  equations.  For  examples,  error  analyses  of  time-delay 
measurements  corrupted  by  noise  are  combined  with  second  stages,  relating 
various  output  errors  to  time-delay  errors,  in  Refs.  5  and  7.  For  tracking 
procedures  in  which  the  relations  between  the  outputs  and  measured  inputs  are 
well  understood  (e.g.,  the  "planar"  problem  in  Ref.  7),  emphasis  is  on  the 
first  stage.  However,  for  a  new  tracking  algorithm  as  considered  here,  some 
basic  issues  associated  mainly  with  the  second  stage  are  essential  to  inves¬ 
tigate.  One  of  these  is  algorithm  robustness,  i.e.,  its  ability  to  yield 
reasonably  accurate  output  information  for  small  imprecision  in  its  inputs. 
Another  related  matter  is  the  determination  of  which  outputs  are  more  sensi¬ 
tive  to  small  or  relatively  large  input  perturbations.  Consequently,  in  one 
portion  of  our  sensitivity  investigation,  we  address  these  issues  by  studying 
propagation  of  non— random  input— measurement  errors,  resulting  from  unspecified 
influences,  on  algorithm  outputs.  We  describe  rather  general  procedures 
for  this  error-propagation  analysis,  and  use  them  in  a  variety  of  tracking 
scenarios.  in  another  portion,  we  determine  relations  between  various 
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quantities  associated  with  random  input  and  output  errors.  We  consider  also 
the  particular  influence  of  noise  on  input  measurements,  to  determine  bounds 
on  tracker  performance  in  terms  of  signal-to-noise  ratio. 

In  Sec.  I,  we  develop  a  nonlinear  system  of  equations  for  the  horizontal¬ 
tracking  algorithm,  relating  the  aforementioned  descriptors  to  parameters 
which  describe  source  position,  motion,  and  frequency.  Section  II  addresses 
the  sensitivity  of  solutions  to  the  algorithm  equations  to  input-measurement 
errors.  An  essential  feature  of  the  methods  discussed  is  a  development  of 
linear  approximations  expressing  algorithm  output  errors  in  terms  of  various 
input  variations.  In  Sec.  Ill,  the  robustness  of  the  algorithm  is  demon¬ 
strated  for  the  propagation  of  input  errors  in  various  source-receiver  scen¬ 
arios.  Some  of  these  scenarios  consist  of  special  cases,  in  which  source  and 
receiver  configurations  are  restricted  in  some  manner.  The  effect  of  random 
input  variations  for  these  special  scenarios  is  discussed  in  Sec.  IV.  The 
influence  of  varying  degrees  of  descriptor -variation  correlation  on  correla¬ 
tion  between  output  errors  is  studied.  We  also  develop  Cramer-Rao  lower 
bounds  for.  the  variances  of  estimates  of  certain  outputs  associated  with  input 
errors  influenced  by  Gaussian  white  noise.  Finally,  major  results  are 
summarized  in  Sec.  V. 

I .  HORI ZONTAL-TRACK ING  ALGORITHM 

We  present  here  a  horizontal-tracking  algorithm  for  the  passive  deter¬ 
mination  of  position,  motion,  and  frequency  of  a  cw  point  source  S,  based  on 
measurements  taken  at  a  linear  receiving  array  R.  The  source  (or  receiver)  is 
submerged  at  constant  depth  Ds  (or  Dr) ,  and  moves  in  a  straight  line  with 
constant  speed  vs  (or  vr) .  The  tracking  algorithm  could  be  adapted  to  more 
general  source  and  receiver  motions,  with  non-linear  paths  and  non-constant 
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speeds,  by  piecewise-linear  approximation.  Source  and  receiver  depths  are 
less  than  about  400  m,  and  the  horizontal  distance  between  S  and  R  is  less 
than  about  10  km.  The  source  is  assumed  to  emit  an  omnidirectional  signal 
which  possesses  a  dominant  single  frequency  fs  (in  Hz) .  The  sound  speed  c  in 
the  ocean  region  between  the  surface  and  max(Ds,Dr)  is  taken  to  be  constant, 
and  conditions  for  the  validity  of  ray  theory  are  assumed.  Thus,  the  source 
signal  travels  in  linear  ray  paths.  For  the  S-R  geometries  considered  here, 
we  regard  the  received  acoustic  field  to  be  dominated  by  the  direct  and 
surface-reflected  rays,  denoted  by  "one"  and  "two",  respectively. 

When  both  rays  arrive  at  R  at  some  time  tn  (in  s) ,  the  direct  ray  is 
emitted  from  S  at  time  t1>n,  while  the  surface-reflected  ray  emanates  at  a 
slightly  different  time,  t2>n.  Figure  1(a)  shows  a  projection  into  a  vertical 
plane  at  time  t^.  Vertical  projections  of  the  one  and  two  rays  are  shown 
here,  along  with  projected  directions  of  S  and  R  (solid  arrows).  Figure  1(b) 
illustrates  a  projection  into  a  horizontal  plane  at  times  t^  and  t2.  The 
horizontal  separation  (range)  between  S  and  R  at  time  tn  is  denoted  by  Rn. 
Angle  3rn  (or  3sn)  is  measured  from  the  horizontal  projection  of  the  path  of 
R  (or  S)  to  the  "line  of  sight"  between  S  and  R  at  tn.  ..Both  3rn  and  &,n  are 
measured  positive  clockwise,  and  are  taken  in  the  interval  -it  <  3^,  <  tt 

[for  example,  in  Fig.  1(b),  3rn  >  0  and  3sn  <  o]. 

Source  location,  motion,  and  frequency  are  determined  at  time  tn  by 
values  of  the  six  parameters  Ds,  Rn,  3rn,  3sn,  vs,  and  fs.  We  note  that 
parameters  R,  3r,  and  3S  vary  with  time,  while  Ds,  vs,  and  fs  have  been  taken 
to  be  time-invariant.  In  order  to  find  these  parameter  values,  we  use  certain 
descriptors  measured  at  R  at,  or  close  to,  two  distinct  times  tn,  n  =  1,2. 
These  descriptors  include  both  directional  information  from  the  linear  array 
and  particular  acoustic  information,  consisting  of  the  period  of  multipath 
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phase  variations  and  the  Doppler-shifted  received  frequency  for  ray  one.  This 
formulation  of  the  tracking  algorithm,  as  well  as  the  sensitivity  analysis  in 
subsequent  sections,  is  extendable  to  the  case  of  more  than  two  measurement 
times.  Utilization  of  more  measurements  would  be  expected  to  yield  improved 
source  localization  and  tracking  estimates.  However,  we  do  not  carry  out  this 
extension  here. 

We  now  develop  equations  at  time  tn  which  relate  our  selected  descriptors 

to  the  six  source  parameters  mentioned  above.  We  shall  present  later  in  this 

section  relations  which  enable  various  parameters  at  time  t2  to  be  expressed 

in  terms  of  those  at  time  t±.  The  directional  information  obtained  from  the 

linear  horizontal  receiving  array  is  the  signal-arrival  angle  for  the  peak 

array  output,  which  corresponds  to  the  direct  ray.  Thus,  the  descriptor  is 

measured  from  the  horizontal  receiver  trajectory  [the  solid  arrow  in  Fig.  1(c)] 

to  the  direct  ray  arrival,  with  0  <  ^  <  it,  and  ^  =  ir/2  being  broadside.  The 

other  line  in  Fig.  1(c)  represents  the  horizontal  projection  of  the  direct  ray 

at  time  tn.  It  is  close  to,  but  not  exactly  the  same  as,  the  line  of  sight  in 

Fig.  1(b),  because  of  the  S  and  R  motion  between  times  t]^n  and  tn.  The  angle 

0n  in  Fig.  1(c)  is  the  vertical  angle  at  time  tn  at  R  of  the  direct  ray,  and 

is  measured  positive  clockwise  from  this  line  to  the  direct  ray  (note  that 

may  assume  positive  or  negative  values,  —  u/2  <  8n  <  tt/2,  depending  on  the 

relative  depths  of  S  and  R)  .  In  Fig.  1(c),  the  angle  0f  at  time  t  is 

rn  n 

measured  in  a  horizontal  plane  from  the  receiver  path  to  the  horizontal  pro¬ 
jection  of  the  direct  ray.  It  can  be  shown  that  0*  equals  the  previously 

rn 

defined  angle  0rn  plus  a  small  correction.  By  approximating  0*  by  0  ,  it 

rn  rn 

follows  that 


cos  =  cos  9n  cos  0rn 


(1) 
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The  quantity  cos  0n  in  Eq.  (1)  can  be  expressed  in  terms  of  the  unknowns  d3 
and  Rn,  so  that 

cos  =  {l  +  [ (Dr”Ds)/Rn]2}  1/2  cos  grn  ,  (2) 

which  is  the  desired  relation  at  time  tn  between  the  descriptor  an^  three 
unknowns . 

We  turn  next  to  the  phase-roll  period,  an  acoustic  descriptor  involving 
multipath  information-  It  has  been  shown1!  for  certain  moving-source/fixed- 
receiver  configurations  that  overall  linear  multipath  phase,  associated  with 
two  dominant  bottom-reflected  rays,  is  interrupted  by  rapid,  regularly-spaced 
changes  in  phase  rate.  These  nearly-periodic  changes  were  shown  to  occur  over 
time  intervals  that  are  from  several  to  tens  of  seconds  long.  If  S  emits  a 
signal  proportional  to  sin(27rfst),  where  t  is  time  in  s,  the  received  signal 
at  R  along  the  direct  (m  =  1)  and  surface-reflected  (m  =  2)  paths  can  be 
expressed  as 

=  Am  sin2Tr[fst  -  <J>m(t)],  m  =  1,2  .  (3a) 

The  relative  amplitudes  Am  can  be  taken  to  be  constant  for  time  intervals  of 
interest.  Although  wm  is  not  a  sinusoidal  function  of  time,  Eq.  (3a) 
expresses  each  arrival  in  the  form  of  one  with  fixed  frequency  f  and  time- 
dependent  phase  <j>m(t)  (in  cycles)  : 

"  ^s  Tmlt''  ”  (m~l)/2,  m  =  1,2  .  (3b) 

In  Eq.  (3b)  ,  xm(t)  is  the  travel  time  of  the  signal  along  the  mth  ray,  and  the 
second  term  accounts  for  a  surface  phase  shift  of  the  m  =  2  ray.  The  total 
acoustic  field  at  R  at  time  t  is 


W(t)  =  A ( t )  sin2ir[fst  -  <f(t)]  , 


(4a) 
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where  the  multipath  relative  amplitude  A(t)  is 

A  ( t)  =  A]_  { 1  +  k2  +  2k  cos  [  A<f>  (t)  j}1/2  ,  (4b) 

and  the  total-field  phase  $  (t)  is 

<fr(t)  =  4>i  ( t )  -  tan”1(k  sin  [  A<j)  (t)  ]  { 1  +  k  cos  [  A<f>(t)  ]  )  .  (4c) 

In  Eqs.  (4b)  and  (4c)  ,  tan“l  is  the  principal  value  of  the  inverse-tangent 
function,  and 


k  =  A2/A1,  A4>  (t)  =  ^(t)  -  4>2(t)  • 


(4d) 


By  differentiating  Eq.  (4c)  with  respect  to  t  and  assuming  (t  -  tn)  to  be 
sufficiently  small  so  that  4>]_(t),  <j>2(t),  and  A<j>(t)  may  be  expanded  to  linear 
terms  in  (t  -  tn) /  we  obtain,  for  the  multipath  phase  rate  (t)  , 


*  (t) 


1  /  n 


-  f  ( x 
s  ^  1  ,n 


(5a) 


where 


G  =  k(k  +  cos  A<j>)(l  +  k2  +  2k  cos  , 


(5b) 


and  x '  =  ( d/dt ) t  (t) | 

m, n  v  J  m  ^=t 


From  Eqs •  (5a)  and  (5b),  $  (t)  repeats  when  A<|> 


changes  by  one  cycle.  Therefore,  the  phase-roll  period  P^  in  s/cycle  of  $  is 


P  =  (f  |  T  -  T  |  ) 
n  ^  s'  1 , n  2 ,  n 1  ; 


-1 


(6) 


where  Eq.  (6)  follows  from  the  linear  [in  (t-t  )  ]  approximation  to  Ab.  By 

2  -1 

using  trigonometry  and  expansions  of  order  M  Rc  and  smaller,  where  Mach 

number  M  is  maxfv  ,v  1/c,  accurate  approximations  for  x  and  x*  mav  be 
s  r  l,n  2,n  *  - 

found.  We  employ  these  approximations  in  Eq.  (6)  to  obtain  an  equation 
relating  p^  at  time  t^  to  the  six  S  parameters  described  previously: 


-9- 


(7a) 


where 


(7b) 


Equations  (6)  and  (7)  suggest  a  connection  of  the  phase-roll  period  to  a 


dif ferential-Doppler  effect,4  but  we  shall  not  discuss  the  similarities  and 


differences  between  these  phenomena  here. 

The  final  descriptor  of  interest  is  the  received  frequency  of  the  direct 
ray  arrival  at  time  tn.  It  is  obtained  by  examining  the  power  spectrum  at  R 
resulting  from  the  emitted  line  spectrum  (at  frequency  fs) .  The  finite-time 
power  spectrum  H(f)  of  the  multipath  arrival  W(t)  is  found  as  a  function  of 
frequency  f  in  the  usual  manner,  with  a  time  average  over  an  interval  of 
length  T,  where  T  is  up  to  several  tens  of  seconds,  about  t  =  tn.  It  can  be 
shown  that  H(f)  has  two  distinct  peaks,  corresponding  to  the  m  =  1  and  2  ray 
arrivals.  Moreover,  H  assumes  the  larger  of  its  t'-'o  peaks  at  frequency 
f  =  frn  for  m  =  1,  since  its  amplitude  is  larger  than  that  for  m  =  2.  The 
quantity  frn  is  given  by 


f 


(8) 


rn 


t  in  Eq.  (8)  is  a  Doppler-type  shift  in  frequency  of  the  direct 
o  x  /  n  * 


and  f  r. 


arrival  resulting  from  relative  source-receiver  motion.  Using  the  approxima¬ 


tion  for  x  discussed  before  Eq.  (7a) ,  we  obtain 


f  =  f  1  +  c 


r  n  s 


(9) 


where  b^/n  is  defined  by  Eq.  (7b) . 


-10- 


We  have  developed  the  three  Eqs.  (2),  (la),  and  (9) ,  which  relate  three 


descriptors  Tpn,  Pn,  and  frn  at  time  tn  to  six  unknowns  Ds/  Rn,  n ,  vs, 


and  fs.  To  supplement  these  three  equations  at  one  time  t^,  we  propose  that 
they  be  employed  at  an  additional  measurement  time  t2.  For  this  purpose,  we 
make  the  reasonable  assumption  that  S  and  R  motion  is  linear  over  time  inter¬ 
vals  between  one  and  several  minutes*  As  discussed  previously,  measurements 
of  the  descriptor  values  have  been  assumed  to  require  shorter  intervals,  of  up 
to  tens  of  seconds. 

The  linearity  of  the  S  and  R  tracks  impose  relations  among  geometric 
quantities  Rn,  Ssn,  and  Brn  at  successive  times.  These  relations  produce  a 
system  of  six  equations  in  six  unknowns  (with  n  =  1)  .  The  solution  to  this 
system,  combined  with  these  relations,  permits  specification  of  5  parameters 
at  t2  as  well.  We  let  At2  =  t2  ~  t]_,  Y2  be  the  length  of  the  line  joining  the 
horizontal  projections  of  R  at  time  and  of  5  at  time  t2  [see  Fig.  1(b)], 
and  n 2  be  the  angle  shown  in  Fig.  1(b).  It  can  be  shown  that  the  following 
equations,  relating  R2 ,  Ss2 9  and  $r2  to  various  n  =  1  unknowns,  are  valid: 


where 


(10b) 


cos  =  [*2  +  R1  ’  (VsAt2)2]/2RlY2 


(10c) 


and 


sin  n. 


'2  =  (VsAt2  Sin  0sJ/Y2  ? 


(10d) 


(11a) 
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and 


r  2 


"  Bs2  +  Pr1 


-  0 


Si 


(12) 


We  omit  tedious  details  of  the  derivation  of  Eqs .  ( 1 Oa ) — ( 1 2 )  for  various  cases 
of  gr1  and  3 s -j  •  These  equations  do  remain  valid  if  the  subscript  2  is 
replaced  by  n,  which  enables  extension  of  our  tracking  procedure  to  more  than 
two  measurement  times. 

In  summary,  the  horizontal  tracking  algorithm  involves  the  solution  of 
the  system  of  six  nonlinear  Eqs.  (2),  (7a),  and  (9)  with  n  =  1,2,  in  terms  of 
the  six  unknowns  R1 ,  3s1 ,  Br1 ,  Ds,  vs,  and  f s .  Equations  (10)  through  (12) 
are  used  to  express  R $$2*  and  $r2  in  terms  of  these  unknowns.  The  solution 
to  this  system  for  typical  input  values  can  be  obtained  numerically  using,  for 
example,  a  nonlinear  least  squares  procedure. we  do  not  provide  an  inves¬ 
tigation  of  existence  and  uniqueness  of  solutions  to  this  system.  This  is 
because  our  principal  interest  here  is  in  algorithm  sensitivity,  and  also 
because  implementation  o^  such  an  algorithm  could  employ  usefully  measurements 
at  additional  times  and  equations  for  S/R  paths  consisting  of  multiple  linear 
segments.  Instead,  we  provide  in  Sec.  Ill  some  typical  examples  of  solutions 
to  the  system  for  various  S/R  scenarios.  After  solution  at  time  t^  is  found, 
Eqs.  (10) -(12)  are  again  employed  to  determine  source  location  parameters  at 
the  later  time  t2 . 


II.  ALGORITHM  SENSITIVITY  ANALYSIS 


In  this  section,  we  present  methods  for  analyzing  the  propagation  of 
measurement  errors  in  horizontal-tracking  algorithm  inputs  on  output  values 
(i.e.,  source  location,  motion,  and  frequency  parameters).  The  inputs  con¬ 
sist  of  the  descriptors  ^n,  Pn,  and  frn,  n  =  1,2,  and  the  sound-speed 
parameter  c.  The  Introduction  described  various  influences  on  measurement 
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errors  that  can  arise  in  tracking  situations.  The  procedures  here  will  be 
used  to  study  output  errors  without  specifying  particular  input-error  causes. 
Input  errors  can  either  be  deterministic  or  random  in  nature,  such  as  those 
due  to  neglect  of  environmental  features  or  noise,  respectively.  Our  pro¬ 
cedures  are  developed  first  for  non-random  types  of  measurement  errors,  so 
that  they  can  be  applied  in  Sec.  Ill  to  study  algorithm  robustness  and  rela¬ 
tive  sensitivities  of  output  errors  to  input  variations.  Later  in  this 
section,  we  present  methods  to  be  used  in  Sec.  IV  for  random  input  variations. 

We  concentrate  in  this  paper  on  examination  of  particular  cases  when  only 
certain  sets  of  input  quantities  are  perturbed:  specifically,  (1)  phase-roll 
periods  and  received  frequencies,  (2)  signal-arrival  angles,  and  (3)  sound 
speed.  The  first  set  consists  of  quantities  that  are  inherently  related  to 
travel-time  variations  due  to  source-receiver  motion  (so  it  is  possible  to 
interpret  errors  in  them  in  terms  of  travel-time  errors).  The  second  set  is 
comprised  of  measurements  from  the  linear  array,  and  errors  in  them  could 
arise  from  various  types  of  array  biases.  The  third  case  is  a  source  of 
ocean-environmental  error  in  this  model.  We  emphasize  that  our  procedures  are 
sufficiently  general  to  apply  to  other  combinations  of  perturbed  inputs* 
including  all  the  descriptors  being  perturbed  simultaneously.  However*  focus¬ 
ing  on  the  particular  cases  clarifies  the  sensitivity  characteristics  of  the 
algorithm. 

We  denote  by  an  overbar,  and  reference  as  "unperturbed”,  all  input  quan¬ 
tities  with  values  obtained  from  a  set  of  measurements,  and,  also,  all  corres¬ 
ponding  output  values  determined  from  the  tracking  algorithm.  All  these 
unperturbed  values  play  the  role  of  an  algorithm  solution  whose  sensitivity  to 
input  perturbations  is  to  be  studied.  Similarly,  "perturbed"  inputs  and  out¬ 
puts  are  unbarred.  Further,  we  use  A  to  indicate  input  perturbations  (for 
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example,  AP]_  =  P]_  -  P]_  is  a  perturbation  in  phase-roll  period  at  time  tj_)  ,  as 
well  as  output  errors  (for  example,  AR2  =  ^2  “  R2  represents  the  error  in 
range  at  time  t2) •  We  also  designate  the  output  error  vector  h^,  n  =  1,2,  by 

^n  =  (ADSr  ARn,  ^^rnf  ^^sn'  ^vs,  Afs)  •  (13) 

Finally,  for  source  depth,  range,  and  source  speed,  the  relative  output  error 
is  defined  to  be  the  output  error  divided  by  the  unperturbed  value,  and  is 
denoted  by  a  x  followed  by  the  appropriate  parameter  in  parentheses  [for 
example,  x(Ds)  =  ADS/DS  is  the  relative  error  in  source  depth]. 

We  now  proceed  to  derive  linear  approximations  for  output  errors  in  terms 
of  input  perturbations.  These  approximations  are  numerically  efficient  in 
that  the  nonlinear  system  of  tracking  equations  only  has  to  be  solved  once 
with  unperturbed  inputs.  They  are  useful  also  for  obtaining  tractable  analy¬ 
tical  formulas  in  the  subsequent  study  of  both  non-random  and  random  input 
variations.  Finally,  although  they  are  derived  only  for  small  input  varia¬ 
tions,  numerical  examples  in  Sec.  Ill  show  them  to  be  accurate  for  rather 
large  input  variations,  as  ^ell.  We  first  develop  linear  approximations  for 
nondimensional  output  errors  in  terms  of  ncndimensional  input  perturbations. 
These  approximations  are  easily  converted  to  dimensional  ones,  are  themselves 
useful  numerically,  and  are  directly  applicable  later  to  comparisons  of 
variances  and  covariances  for  random  output  errors  and  input  variations.  Non- 
dimensional  quantities  using  unperturbed  input  values  are  defined  as  follows; 

Pn  =  Rn/Rl '  ^n  =  ^rn/^rl'  n  =  1*2,  (14a) 

and 

C  =  c/c,  e  =  (Pifrl)”1'  ~  =  (Pi.fr l”r )/°  • 


(14b) 
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Dimensionless  unknowns  Xjn,  for  1  <  j  <  6  and  n  =  1,2,  are  defined  to  be 
xln  =  C^r“Ds)/Rn'  x2n  =  (Dr+Ds)/Rn'  x3n  =  &rn'  x4n  =  ^sn' 


x5n  =  C^lfrlvs )/c'  x6n  =  fs/frl 


Note  that  Ds  and  Rn  can  be  recovered  from  X]_n  and  X2n  by 


Ds  “  Dr  (x2n~xln)/(xln+x2n)  •  Rn  ~  ^Dr/ (x]_n+X2n ) 


(14c) 


(14d) 


A  nondimensional  version  of  the  tracking  algorithm,  Eqs.  (2) ,  (7a) ,  and  (9) , 

is  given  for  n  =  1,2  by 


(i 


+  x 


2  rl/2 

lnJ 


cos  x^  -  cos  =  0  , 
3n  Yn 


(15a) 


Cx6ih(l+xin)’V^  -  (1+x2n)”VZKx5nCOS  x4n+v  cos  x3rJI 


2  >,-1/2 


2  >,-1/2- 


i-l 


-  Pn  "  0 


(15b) 


and 


••  “1  2  1/2 
x.  [  1  +  e  C  fl+x-,  1  /  fx_  cos  x„  +  v  cos  x.  1 1  -  f  =  0. 

fin  L  ^  n-'  '•fin  4n  '  J  n 


6n 


Nondimensional  output  vectors  xn  are 


(15c) 


=  fxn  ,x„  ,x_  ,x.  ,x^  ,x  _  1  ,  n  =  1,2  , 
^  In  2n  3n  4n  5n  6n } 


(16) 


and  from  Eqs.  (10)-(12),  formulas  for  the  components  of  xr>  (not  written  here) 
can  be  obtained  in  terms  of  those  of  x^ .  Finally,  the  nondimensional  input- 
perturbation  vector  5,  which  can  represent  any  of  the  three  sets  of  perturba¬ 
tions  described  previously,  is  defined  by 


6  =  (  *  <$2 '  ^3  '  ^4  '  ^5 '  ^6 )  • 


(17a) 
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For  each  set,  component  6j_ ,  i  =  1,6,  corresponds  to  the  perturbation  in  the 
ith  equation  of  Eqs.  (15),  written  with  the  n  =  1  equations  preceding  the 
n  =  2  equations.  For  example,  with  nondimensional  perturbations  in  phase-roll 
period  Apn  and  received  frequency  Afn,  n  =  1,2,  we  have 

6  =  ( 0,Ap1 ,Af 1 ,0,Ap2,Af2)  .  (17b) 

Linear  approximations  for  the  nondimensional  error  vectors  Axn(n=1,2), 
applicable  to  all  input-variation  sets,  are  now  written  as  follows.  We 
express  Eqs.  (15)  as 

=  o,  i  =  1,6  ,  (1?c) 

and  let  J  be  the  6x6  Jacobian  matrix  for  Eqs.  (17c).  Then  entry  (i,j)  of  J 
is 

Jij  =  ^i/^xjl  /  (17d) 

where  the  partial  derivatives  are  evaluated  at  (x1#5)  =  (x>j,o),  and  0  is  the 
zero  vector.  We  do  not  list  the  entries  of  J  in  this  paper.  We  let  E  be  the 
6x6  matrix,  the  kth  column  of  which  consists  of  partial  derivatives  of  the 
Fj_  with  respect  to  <5k  evaluated  at  (^,5)  =  (x^,o)*  For  all  three  sets  of 
input  perturbations,  the  matrix  E  is  diagonal.  It  can  be  shown  by  generaliz¬ 
ing  a  scalar  result21  that  a  linear  approximation  to  Ax-j  ,  valid  to  order  6,  is 
• 

ASl  =  ,  (18a) 

where 

B1  =  -J_1E  (18b) 

and  J-1  is  the  inverse  of  the  Jacobian  matrix.  We  note  that  Eqs.  (18)  may  be 
solved  numerically  in  an  efficient  manner  by  a  factorization  scheme20  applied 
to 
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J(A5l)  =  “  E§  •  (19) 

Numerical  experience  with  Eq.  (19)  indicates  that  it  generally  produces  more 
accurate  solutions  than  a  corresponding  dimensional  linear  system  derived  from 
Eqs •  ( 2) ,  (7a),  and  (9).  By  using  Eq.  (18a),  the  expressions  mentioned  prev¬ 
iously  for  x2  in  terms  of  x^,  and  a  linear  expansion  of  Ax£,  we  obtain 
similarly  that 

Ax2  =  B2  6  ,  (20a) 

where 

b2  =  K  Bi  .  (20b) 

The  6x6  matrix  K  consists  of  partial  derivatives  of  the  components  of  Ax^ 
evaluated  at  Ax^  =  0,  and  B]_  is  given  by  Eq.  (18b).  From  Eqs.  (14c),  (14d)  , 

(18a) ,  and  (20a) ,  approximations  for  the  dimensional  error  vectors  h^  of 
Fq.  (13)  are  readily  obtained: 

”  L^Br(xln  +  x2n)  ^(xln^x2n  ~  x2n^xln]'  "  ^^r  C^ln  +  x2n  )  ^  ^ 

(Axln  +  Ax2n),  Ax3n,  Ax4n,(c/P1Erl)Ax5n/  frlAx6n],  n  =  1,2  .  (21) 

Approximations  for  various  relative  output  errors  x  follow  immediately  from 
Eq.  (21). 

We  turn  briefly  to  formulas  useful  in  analyzing  output  errors  that  are 
associated  with  random  input  variations.  All  input  perturbations  and  output 
errors  are  modeled  as  random  variables,  but  without  specification  of  their 
probability  distributions.  These  random  variables  possess  mean  values  of 
zero.  Equations  (18a)  and  (20a)  imply  that  the  nondimensional  output-error 


random  variables  satisfy 
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A?n  =  Sn  ^  n  =  1,2  ,  (22) 

where  the  input  error  S  is  a  vector  random  variable,  but  the  entries  of  Bn 
are  non-random.  Entry  (i,j)  of  the  6x6  error  covariance  matrix^  V 
associated  with  5  is 

vij  E  (6i6k)  »  (23a) 

where  is  expectation.  A  similar  equation  holds  for  the  error  covariance 
matrix  Vn  associated  with  Axn.  It  can  be  shown^2  that  V  and  Vn  satisfy 

Vn  «  ^  V  Bh  ,  (23b) 

where  -  denotes  matrix  transpose.  Equation  (23b)  is  the  desired  relationship, 
expressing  output-error  variances  and  covariances  in  terms  of  those  of  input 
perturbations • 

III.  ERROR  PROPAGATION  ANALYSIS 

We  analyze  the  sensitivity  of  the  horizontal-tracking  algorithm  outputs 
to  deterministic  input-measurement  errors  in  a  variety  of  source-receiver 
configurations.  Typical  scenarios  are  displayed  in  top -views  in  Fig.  2. 
Scenario  (A)  shows  S  and  R  on  headings  which  are  perpendicular  at  time  t^ ,  so 
that  the  range  R  between  them  shortens  at  a  later  time  t2.  Scenario  (B)  is  an 
example  of  diverging  range,  while  S  and  R  converge  at  varying  rates  in  (C), 
(D),  and  (E) .  Unperturbed  values  for  algorithm  parameters,  descriptors,  and 
track  outputs  are  given  in  Table  I*  In  subsection  A,  we  consider  scenarios 
(A)-(C);  in  subsection  B,  we  discuss  scenarios  (D)  and  (E),  in  which  motion  is 
confined  to  a  vertical  plane,  along  with  a  particular  version  of  scenario  (A). 
We  remark  that  we  considered  numerous  other  cases  as  well,  which  are  not 


reported  here. 
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A.  Example  scenarios 

For  each  of  scenarios  (A)-(C),  track  outputs  were  computed  using  the 
numerical  procedure  described  in  Sec.  I  for  the  full  nonlinear  system  of 
algorithm  Eqs.  (2) ,  (7a) ,  and  (9) .  The  principal  result  of  this  subsection  is 
that,  for  all  scenarios,  small  perturbations  of  all  algorithm  inputs  produce 
small  variations  in  output  estimates.  For  example,  we  display  in  Fig.  3(a) 

[or  Fig.  3(b)]  graphs  of  relative  percentage  errors  in  depth,  range,  and  speed 
at  time  t2  versus  phase-roll  period  perturbation  Ap(  =  AP]_  =  AP2  ),  for  scenario 
(A)  [or  (B )  ]  .  The  value  of  frequency  perturbation  Afr(  =  Afri  =  Afr2  )  is 
0.1  Hz,  while  AP  varies  from  0  to  3  s/cycle.  These  near ly-linear  curves  in 
Figs.  3(a)  and  3(b)  show  that  output  errors  are  small,  not  only  for  small 
values  of  AP,  but  also  for  relatively  large  values.  For  example,  Fig.  3(a) 
indicates  that  relative  output  errors  are  less  than  8%  for  AP  values  corres¬ 
ponding  to  relative  period  errors  up  to  about  20%.  The  speed  error  (or  depth 
error)  curve  is  the  largest  in  Fig.  3(a)  [or  Fig.  3(b)  ],  which  shows  that  the 
relative  magnitudes  of  output  errors  are  dependent  upon  scenario.  An  example 
of  an  approximate  error  curve,  computed  using  Eq.  (21),  is  displayed  in 
Fig.  3(b)  as  the  upper  long-dashed  curve  for  range  error*  It  shows  the 
accuracy  of  the  linear  approximations  to  output  errors,  for  both  small  and 
relatively  large  values  of  AP.  Also,  for  scenarios  (A)-(C),  output  errors 
were  shown  to  be  much  less  sensitive,  to  changes  in  received  frequencies  than 
to  phase-roll  period  perturbations. 

All  output  predictions  are  relatively  insensitive  to  small  input  varia¬ 
tions.  However,  some  outputs  are  more  sensitive  to  larger  variations  in  some 
inputs  than  in  other  inputs.  For  example,  we  show  in  Fig.  4  relative  depth, 
range,  and  speed  percentage  errors  at  time  t2  in  scenario  (C) .  Figure  4(a) 
illustrates  these  output  errors,  with  an  ordinate  scale  up  to  25%,  versus 
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AP  (with  Afr  =  0.1  Hz).  Figure  4(b)  shows  the  same  errors,  with  ordinate  up 
to  5%,  versus  Aip  (=  A^]_  -  A^)  /  where  Aip  represents  the  measurement  error  of 
the  signal-arrival  angle.  Results  in  Fig.  4  are  typical  for  other  scenarios, 
as  well. 

We  turn  now  to  the  influence  of  various  input  perturbations  on  errors  in 
bearing  and  course.  In  Fig.  5(a)  [or  Fig.  5(b)],  graphs  of  bearing  -  and 
course  -  angle  errors  at  time  t2  versus  A  if;  are  displayed  for  scenario  (B)  [or 
scenario  (C) ] .  Results  show  that  bearing  (A0r)  and  course  (A8g)  angle  errors 
are  comparable  to  signal-arrival  angle  perturbations.  We  found  this  conclu¬ 
sion  to  be  widely  applicable.  It  follows  that  the  main  influence  of  Aip  or 
A8r  occurs  in  algorithm  Eq.  (2).  These  errors  in  Sr  influence,  in  turn,  the 
errors  in  3S  in  Eqs.  (7a)  and  (9),  where  both  these  quantities  are  present. 

In  Fig.  5(a)  [or  Fig.  5(b)]  values  of  A8r2  are  always  greater  (or  less)  than 
A6S2*  which  shows  that  relative  error  magnitude  varies  with  scenario.  In  con¬ 
trast  to  the  comparable  influence  of  Aip  on  both  A8r  and  Ags/  there  is  a 
significant  difference  in  the  effect  of  AP  on  these  quantities.  Small  varia¬ 
tions  in  phase-roll  period  yield  small  changes  in  A0r  and  Agg,  but  larger 
values  of  AP  produce  significantly  larger  values  in  Ags__than  in  Agj-  in  all 
scenarios.  For  example,  with  scenario  (C)  at  time  t2  and  values  of  AP  = 

2  s/cycle  (and  Afr  =  0.1  Hz),  1 ABS2 1  =  3.0°,  while  |  A3r2 1  =  0.05°.  Since  A\|/ 
is  zero  here,  the  relation  between  8r  and  ip  in  Eq.  (2)  evidently  serves  to 
diminish  the  influence  of  AP  on  A8r  in  Eq.  (7a). 

The  input  perturbation  which  least  affects  output  errors  is  As.  For  Az 
values  up  to  10  ms-^  in  scenarios  (A)-(C),  output  errors  are  nearly  zero. 

B.  Special  tracking  scenarios 

We  discuss  here  the  sensitivity  of  output  errors  to  input  perturbations 
in  particular  types  of  tracking  scenarios.  Scenarios  (D)  and  (E)  in  Fig.  2 
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have  motion  confined  to  a  vertical  plane,  and  we  consider  a  variant  of 
scenario  (A)  with  perpendicular  tracks  at  time  t^ .  These  scenarios  are 
"special”  in  that  a  priori  knowledge  is  assumed  for  bearing  angles  ^n,  n  = 
1,2,  course  angles  3sn,  n  =  1,2,  and  source  frequency  fs.  For  subsequent 
convenience,  we  also  take  receiver  depth  less  than  source  depth.  The  moti¬ 
vation  for  investigating  special  tracking  scenarios  is  to  obtain  tractable 
analytical  expressions  for  output  errors  in  terms  of  input  perturbations. 
Furthermore,  useful  information  may  be  obtained  from  these  for  both  determin¬ 
istic  input  variations  in  this  subsection  and  for  random  variations  in  Sec. 

IV.  Finally,  certain  insights  into  the  robustness  of  the  tracking  algorithm 
in  more  general  scenarios  are  gained  by  understanding  sensitivity  results  in 
special  tracking  scenarios. 

The  a  priori  knowledge  for  special  scenarios  is  useful  to  reduce  the  full 
tracking  system,  Eqs.  (2) ,  (7a) ,  and  (9) ,  to  three  tractable  nonlinear  equa¬ 
tions  at  time  t]_.  In  this  reduced  system,  descriptors  ,  and  frj_ 

determine  three  source  parameters  Ds,  R]_,  and  vs  via  explicit  formulas.  The 
remaining  unknown  R2  is  then  found  from  Eq.  (10a) .  These  formulas  will  not  be 
presented  here.  Instead,  the  accurate  linear  approximations  corresponding  to 
Eq.  (21)  will  be  described.  It  can  be  shown  that  Eqs.  (15)  can  be  reduced  to 
a  nondimensional  system  of  two  equations  in  two  unknowns  X21  and  x5^.  By 
applying  Eq.  (21)  to  this  system,  approximations  for  output  errors  are 
obtained.  For  examp-le,  it  can  be  shown  that,  for  perturbations  in  period 
and  frequency  Afr^,  relative  errors  x(Ds)  '  x(R]_)  /  and  X(vs^  are 


X(DS)  =  [ (Ds-Dr)/Ds]c1(Y1+y2)  , 

(24a) 

X  ( R1 )  =-C1(Yi+Y2) 

(24b) 

and 
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X(vs)  =  (cos  esl/cos  erl)[c(frl-fs)/vs  fs  cos  4*1  ] Y2  ,  (24c) 

where  the  relative  period  and  frequency  variations  are 

Yi  =  APf/Pi  (24d) 

and 

y2  =  Afrl/(frl-fs)  .  ( 24e) 

The  quantity  c^  is 

C1  =  (Rq/2Dr)(cos  3r]/cos  ^i) {02/ [l- (03  cos  *l)2F2}  .  (25a) 

where 

cl  E  [l-Cfs/fri)]/{Pi  frlt^ffs/frl)  -  (l/Pi  fri)]2}  (25b) 

and 

c3  =  [l-(Vfrl)  -  (VPifri)]/[l-(fs/fri)]  .  (250 

Values  of  Ds,  Rj.,  and  vs  in  Eqs.  (24) -(25)  are  found  from  the  reduced  system 
described  at  the  start  of  this  paragraph.  We  conclude  from  Eqs.  (24a)  and 
(24b)  that  under  our  assumption  Dr  <  Ds,  the  magnitudes  of  relative  errors  in 
range  are  greater  than  those  in  source  depth.  Also,  we  note  that  relative 
error  in  source  speed  is  influenced  only  by  the  received  frequency  perturba¬ 
tion,  unlike  relative  errors  in  depth  and  range  at  time  ti  . 

For  scenarios  (A) ,  (D) ,  and  (E) ,  we  found  that  typically  small  input  per¬ 
turbations  resulted  in  small  errors  in  all  outputs.  This  result  is  entirely 
consistent  with  results  in  Sec.  Ill  A  for  more  general  S-R  configurations. 

For  example,  we  display  in  Fig.  6(a)  [or  Fig.  6(b)]  graphs  of  percentage 
relative  errors  in  source  depth,  speed,  and  range  at  times  tj_  and  t-  versus 
AP  =  AP]_  for 'special  scenario  (A)  [or  scenario  (D)  ]  ,  with  Afrl  having  the 
constant  value  0.1  Hz.  We  show  in  Figs.  6  and  7  graphs  obtained  from  Eqs. 
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(24)— (25)  and  from  an  approximation  (not  included  explicitly  here)  for  x(^2^ • 
Numerous  numerical  calculations  indicate  the  good  accuracy  of  these  approxima¬ 
tions.  Figs.  6(a)  and  6(b)  demonstrate  that  output  errors  are  small  for  both 
small  and  relatively  large  values  of  AP.  A  comparison  of  Figs.  3(a)  and  6(a) 
indicates  that  these  percentage  errors  are  somewhat  larger  for  the  restricted 
scenario  (A)  than  for  the  general  case  with  scenario  (A) .  We  observe  also 
various  qualitative  conclusions  of  the  output  error  approximations  in  Figs. 
6(a)  and  6(b)  as,  for  example,  speed  errors  in  both  figures  being  independent 
of  AP  perturbations. 

In  the  previous  subsection,  output  errors  were  affected  much  more  signif¬ 
icantly  by  AP  than  by  Afr  perturbations.  For  the  special  scenarios  here, 
frequency-measurement  errors  can  influence  output  errors  comparably  to  phase- 
roll  period  variations.  This  can  be  seen  from  Eqs.  (24)  ,  since  relative 
period  and  frequency  errors,  and  y 2 /  assume  values  of  comparable  magnitude. 
Also,  Fig.  7(a)  [or  Fig.  7(b)]  shows  the  effect  of  Afr  values  up  to  1  Hz  (with 
AP  =  1  sec)  on  percentage  depth,  range,  and  speed  errors  for  restricted 
scenario  (A)  [  or  scenario  (E) ] .  Similarly,  the  influence  of  the  received 

angle  perturbation  on  percentage  errors  in  depth  and  range  is  considerably 
larger  for  the  special  scenarios  than  in  subsection  A. 

Our  analysis  in  this  subsection  has  corroborated  the  conclusions  in  Sec. 

I IIA  in  that  the  tracking  algorithm  is  generally  robust  for  reasonable  input 
perturbations.  The  investigation  here  was  performed  on  a  subsystem  of  the 
equations  associated  with  the  full  algorithm.  It  also  indicates  that  for 
larger  input  perturbations,  certain  output  errors  are  typically  much  less  with 
the  full  algorithm  than  with  an  algorithm  based  on  the  subsystem.  Therefore, 
the  full  algorithm  is  better  able  to  absorb  larger  input  fluctuations  and 
still  produce  good  predictions.  This  suggests  enhanced  robustness  of  the  full 
algorithm  by  incorporating  measurements  at  more  than  two  times. 
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IV.  SENSITIVITY  TO  RANDOM  VARIATIONS 

We  discuss  relations  between  random  input  measurements  and  corresponding 
output  errors  for  our  tracking  algorithm.  Examples  of  possible  sources  of  ran¬ 
dom  influences  on  input-measurement  errors  are  discussed  in  the  Introduction. 

In  subsection  A,  we  study  random  output  errors  without  specifying  any  particu¬ 
lar  input-error  mechanism;  in  subsection  B,  we  consider  the  specific  effect  of 
input  noise  variations.  In  order  to  present  tractable  formulas,  we  focus  on 
special  tracking  scenarios  as  described  in  Sec.  IIIB.  However,  we  stress  that 
the  methods  used  with  the  special  scenarios  can  be  generalized  to  core  general 
ones.  We  choose  to  analyze  effects  of  phase-roll  period  and  received  frequency 
perturbations,  at  time  tp  and  with  zero  means,  on  relative  errors  in  source 
depth  and  speed.  Random  relative  range  errors  are  not  included,  since  they  can 
always  be  derived  from  relative  source-depth  errors. 

A.  Random  error  correlations 

In  Sec.  II,  it  is  noted  that  the  period  and  received  frequency  are  both 
related  to  travel-time  variations  due  to  S-R  motion.  It  is  reasonable,  then, 
to  expect  that  random  influences  which  affect  errors  in..one  of  these  descrip¬ 
tors  also  produces  errors  in  the  other.  We  wish  to  describe  the  effect  of 
varying  correlation  between  perturbations  in  these  descriptors  on  the  correla¬ 
tion  between  output  errors.  Zero-mean  random  variables  for  nondimensional 
period  and  frequency  variations  are  represented  here  by  yp  and  y2  in  Eqs.  (24d) 
and  (24e)  ,  and  by  x(°s)  and  X(vs)  for  relative  source  depth  and  speed  errors. 

The  input-perturbation  variances  and  covariance  are  a2  =  E{y2),  i  =  1,2,  and 

yi  i 

E(YiY2)'  where  E  is  again  the  expectation  operator.  The  correlation 
coefficient  pp  for  Yq  and  Y2  is 


°r  "  EIYiV/0nV  ' 


(26a) 
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and  a  similar  definition  holds  for  associated  with  x<Ds)  and  X(vs)-  Using 
Eqs.  (24a)  and  (24c)  as  valid  linear  approximations  for  x(Ds^  anc*  X(vs)  in 
terms  of  Y-|  and  y2,  along  with  a  modified  version  of  Eq.  (23b),  the  following 
formula  for  is  derived: 


^  2  1/2 
p  =  ( 1  +  <Pp)/l  1  +  k  +  2kp  )  / 

X 


(26b) 


where 


(26c) 


2 


It  is  important  to  observe  that  p^  from  Eq.  (26b)  is  influenced  solely  by  the 
correlation  coefficient  pp,  and  the  ratio  k  of  standard  deviations  associated 
with  relative  period  and  frequency  perturbations.  That  is,  it  is  independent 
of  all  mean-state  quantities  and  constants  in  Eqs.  (24a)  and  (24c). 

We  display  in  Fig.  8  graphs  of  p^  versus  pp  for  various  values  of  k 
ranging  from  zero  to  infinity.  From  the  k  =  0  and  <  =  0.5  curves,  associated 
with  relatively  small  period-error  standard  deviations,  p^  is  close  to  one. 
This  indicates  that  relative  errors  in  source  depth  and  speed  are  nearly 
perfectly  correlated  (and  nearly  linearly  related^),  regardless  of  the  cor¬ 
relation  between  period  and  frequency  variations o  As  k  increases,  the  curves 
in  Fig.  8  decrease  monotonically .  For  large  values  of  k,  p^,  approaches  the 
diagonal  k  =  00 .  Thus,  when  the  standard  deviation  of  period  perturbations  is 
significantly  greater  than  that  associated  with  frequency,  output-^rror  corre¬ 
lation  is  virtually  identical  with  the  correlation  between  input  variations. 
The  fact  that  all  curves  in  Fig.  8  begin  at  the  points  (pp,p^)  =  (-1 ,  ±1 )  and 
end  at  (1,1)  shows  that  a  perfect  correlation  between  y-j  and  y2  always 
produces  the  same  result  in  x(Ds)  and  x(vs)*  The  particular  case  <  =  1  ,  with 
a  discontinuity  at  pp  =  -1 ,  represents  a  transition  between  a  situation  where 


p^  is  never  zero  (<  <  1 )  and  one  where  p^  is  always  zero  for  some  nonzero 
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value  of  (k  >  1).  Finally,  we  conclude  that  if  pp  =  0  (input  perturbations 
are  uncorrelated),  then  ±  0  (output  errors  possess  some  correlation), 
unless  ic  is  large, 

B.  Tracker  performance  bounds 

Next,  we  regard  the  received  signal,  from  which  phase-roll  period  and 
frequency  measurements  are  obtained,  as  corrupted  by  noise.  The  noise  N(t)  is 

taken  to  be  Gaussian  with  zero  mean,  spatially  incoherent,  and  spectrally 
white  with  power  No  each  one-Hertz  band.  We  express  the  received  signal 
Z ( t )  at  R  at  time  t  as 

Z(t)  =  W(t)  +  N(t)  ,  (27a) 

where  W  is  given  by  Eq.  (4a).  Errors  in  estimates  for  the  period  and  received 
frequency  are  taken  to  possess  Gaussian  probability  densities.  As  before, 
means  of  these  descriptors,  along  with  output  source  depth  and  speed  values, 
are  taken  as  known.  With  t^  =  0  for  convenience,  the  observation  time  T  of 
measurements  extends  from  -  T/2  to  T/2  s,  and  encompasses  tens  of  seconds  (see 
Sec.  I) •  The  time-averaged  power  Sg  of  the  signal  is 23 

T/2  2  22 

S  =  lim  (i/t)  /  W  ( t)dt  =  (a;  +  A^)/2  .  (27b) 

°  -T/2  1  2 

One  appropriate  signal-to-noise  ratio  (SNR)  here  is  then  Sg/Ng.  Cramer-Rao 
lower  bounds24  for  the  variances  of  unbiased  estimators  for  the  track  outputs 
are  desired,  as  functions  of  SNR.  A  two-stage  estimation  procedure  is 
required:  First,  obtain  such  bounds  for  variances  of  descriptors  in  terms  of 

SNR;  subsequently,  express  track  output  variance  bounds  in  terms  of  those  for 
descriptor  variances. 

For  the  first  stage,  the  2-vector  a  is  defined  with  a1  =  p-j  and  ct2  =  f-j 
as  nondimensio‘nal  versions  of  period  and  received  frequency.  Let  ¥(a)  denote 
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the  symmetric  2><2  Fisher  Information  matrix^  associated  with  a.  Applying 
the  approach  for  sinusoidal  signals  in  Ref.  4,  it  can  be  shown  that  j  ( ot) , 
the  ijth  entry  of  ¥(ot),  is  given  by 


T/2 

Yij(a)  ='(1/N0)  / 

-T/2 


Ow/3^)  (  3w/3aj  )dt  . 


(28) 


Entries  of  'F(a)  are  evaluated  from  Eq.  (28)  by  representing  W  as  the  sum  of 

the  two  signals  in  Eq.  (3a) ,  expressing  their  phases  <j)-j(t)  and  <j>2(t)  in  terms 

of  a-j  and  a2,  and  using  the  relatively  large  size  of  T  to  neglect  small  terms. 
Cramer-Rao  lower  bounds  for  the  variances  a^(a^)  and  a^( a2)  of  unbiased 

estimators  of  and  a2  are  the  diagonal  elements  A -|  -j  ( ct )  and  A22(a)  of  th e 

symmetric  2x2  matrix  A(a)  =  [¥(a)]”1.  It  follows  that 

q2(S1)  >  A-|  -j  (a)  =  (P1/2ir)2(12/T3)(k  +  k'1 )  2(  S0/N0) -1  (29a) 

and 


a2(a2)  >  A22(a)  =  0/2*  fr1 )  2(  1  2/T3 )  ( 1  +  k2)  ( S0/N0) -1 


(29b) 


where  k  is  given  in  Eq.  (4d).  We  also  require  later  the  off-diagonal  entry, 


A12(?)  =A21(?)  =  -  (Pl/2lT  f r  1  )  2(  1  2/T3 )  (  1  +  k"2)  (Sq/Nq)-1  .  (29c) 

We  note  that  Eqs .  (29)  depend  on  Sq/Nq,  T,  k,  and  the  dimensional  mean 
values  of  period  and  received  frequency.  It  follows  from  Eqs.  (28a)  and  (28b) 
that  the  Cramer-Rao  lower  bound  associated  with  period  perturbations  a-j  is 
much  greater  than  that  for  frequency  perturbations  ot2.  Also,  for  our  algor¬ 
ithm  involving  multipath  information  at  a  single  receiver,  bounds  decrease  as 
1/T^  for  increasing  T.  This  qualitative  result  agrees  with  that  reported^  for 
differential  Doppler  estimation  for  single-path  signals  at  multiple  receivers. 

Turning  to  the  aforementioned  second  stage,  we  establish  Cramer-Rao  lower 
bounds  associated  with  the  2-vector  z,  where  z-j  =  Ds/Dr  and  z2  =  X5-]  [from 
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Eq.  (14c) ]  are  nondimensional  versions  of  source  depth  and  speed.  From  the 
system  of  two  equations  for  and  f]_  =  a2  described  in  Sec.  Ill  b, 

relations  a  =  a(z)  can  be  found.  The  ijth  entry  Q^j  of  the  2x2  Jacobian 
matrix  Q  for  this  new  system  is 

Qj_j  =  9a^/3  z  j  ,  (30) 

where  the  partial  derivatives  are  evaluated  at  mean  values  of  z  components. 
For  our  situation, 

Qll  =  2x61e1|z2  cos  lsl  +  v  cos  0rl |  (l+eie2) _3/2  [  (l+Ii) / (1-Zi)  3  ]  ,  (31a) 
Ql2  =  x6|cos  gsl|[  (cos  Vcos  lrl)  -  (l+e1e2)_1/2]  f  (31b) 

222  =  e  x6 (cos  &si/cos  3ri)cos  ^  ,  (31c) 

where 

e-L  =  (cos  *grl/cos  t!)  2  -  1,  e2  =  [  (M-z^/fl-z^  ] 2  ,  (31d) 


and  e,  v,  and  xgi  are  defined  in  Eqs.  (14b)  and  (14c) 0  Under  our  model 
assumptions,  it  follows5  that  the  Fisher  Information  matrix  ’f(z)  associated 
with  z  is 

f(z)  =  Q  f(a)Q  .  (32) 


Cramer-Rao  lower  bounds  associated  with  z  are  the  diagonal  elements  of  A(z)  = 

tf(z)]"1,  as  before.  From  Eqs.  (29),  (31),  and  (32)  and  the  definitions  of  zi 
and  z2,  such  bounds  for  the  variances  o2(Ds),  a(vs)  of  unbiased  estimators  of 


Ds  and  vs  are 


»2(Bs)  >02  q-2  [AU(S)  -  2(Q12/Q22)A12(«)  ♦  (Q12/822)2A22(a)] 


(33a) 


-28- 


and 


o2{v  )  >  (c/p  f  )2  Q-2  A  fa)  . 
s  1  rl  22  22 


(33b) 


From  Eq.  (33a)  the  lower  bound  on  source  depth  variances  depends  on 
bounds  for  both  period  and  received  frequency  variances.  From  Eq.  (31b),  the 
lower  bound  for  speed  variances  depends  only  on  the  lower  bound  for  received 
frequency  variance.  This  result  is  analogous  to  one  in  Sec.  Ill  B  for  non- 
random  error  variations.  We  display  in  Fig.  9  graphs  of  the  Cramer-Rao  lower 
bounds  for  depth  and  speed  error  for  special  scenario  (A)  versus  the  dB  SNR, 

10  log]_Q  (Sq/Nq)  .  The  lower  bounds  are  shown  as  percentage  standard  devia¬ 
tions,  normalized  via  the  mean  value  of  each  output.  Measurement  observation 
time  used  for  both  curves  in  Fig.  9  is  T  =  30  s.  The  very  small  values  for 
the  speed  bound,  of  less  than  about  2%,  corresponding  to  a  speed  error  of  less 
than  0.1  ms~-*-,  are  noteworthy.  They  arise  because  of  its  aforementioned 
dependence  on  the  received  frequency  lower  bound.  This  bound,  ^2  (  °0 '  is  very 
small  due  to  the  quantity  (1/2tt  fri)  2  in  Eq.  (29b).  The  depth  bound  in  Fig.  9 
decreases  steadily  from  just  under  30%,  or  90  m,  at  a  SNR  of  0  to  just  under 
10%,  or  30  m,  when  the  SNR  is  10.  This  bound  is  influenced  mainly  by  the 
phase-roll  period  variance  bound,  and  is  consequently  larger  than  that  for  the 
speed  curve.  Results  for  scenarios  (D)  and  (E)  are  very  similar  to  those 
shown  in  Fig.  9.  We  conclude  that  variance  bounds  for  the  special  scenarios 
are  excellent  for  source  speed  and  good  for  source  depth  at  modera-te  to  high 
SNR. 

Our  procedures  could  be  extended  to  analyses  of  the  more  general  scen¬ 
arios.  As  illustrated  in  Sec.  Ill,  output  errors  tended  to  be  comparable  to 
or  less  than  those  associated  with  the  scenarios  considered  here.  Thus,  it  is 
anticipated  that  bounds  on  the  performance  of  the  full  algorithm  should  be 


reasonable. 
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V.  SUMMARY 

This  paper  considers  the  passive  determination  of  bearing,  range,  depth, 
course,  speed,  and  frequency  for  a  moving  cw  acoustic  source  S.  Measurements 
are  taken  from  a  moving  linear  receiving  array  R.  Both  S  and  R  are  submerged 
at  constant  depths  of  less  than  400  m,  separated  by  ranges  less  than  about  10 
km,  and  they  move  on  linear  paths  at  constant  speeds.  The  dominant  acoustic 
signals  are  assumed  to  travel  along  two  upper-ocean  ray  paths,  one  direct  and 
one  surface-reflected,  and  the  sound  speed  is  taken  to  be  constant  there.  The 
tracking  is  accomplished  using  a  new  combination  of  acoustic  quantities  (or 
descriptors) .  These  include  the  signal-arrival  angle  at  R  and  the  Doppler- 
shifted  received  frequency  that  are  associated  with  the  direct  ray,  as  well  as 
the  phase-roll  quasi-period  of  multipath  phase  variations.  The  latter 
involves  regularly  spaced,  brief  changes  in  total-field  phase  rate.  Using 
descriptor  measurements  at  two  distinct  times,  a  system  of  equations  for  the 
horizontal  tracking  algorithm  is  developed,  relating  descriptors  to  six  un¬ 
known  source  parameters.  Our  algorithm  formulation  is  sufficiently  general  so 
that  it  can  be  extended  to  more  than  two  measurement  times  for  additional 
improvement  in  source-parameter  estimates «  Solution  examples  via  a  least- 
squares  numerical  procedure  are  provided  for  various  scenarios. 

We  then  investigate  the  sensitivity  of  output  values  to  variations  in 
inputs,  which  include  the  descriptors  and  the  sound  speedo  Such  variations 
might  result  from  measurement  errors  associated  with  noise,  the  receiving 
array,  environmental  fluctuations,  and  other  factors.  Procedures  are  pre¬ 
sented  for  analyzing  the  effects  of  both  deterministic  and  random  input 
perturbations  on  output  errors.  The  methods  rely  on  efficiently  octaining 
linear  approximations  for  output  errors  in  terms  of  input  variations. 


These 
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approximations  are  also  computationally  efficient,  and  results  using  them  are 
shown  to  be  valid  for  small,  as  well  as  somewhat  larger,  variations*  Further, 
they  are  useful  to  express  variances  and  covariances  of  localization  errors  in 
terms  of  those  for  random  input  variations* 

We  focus  on  three  typical  cases  of  perturbed  inputs,  associated  with 
phase-roll  periods  and  received  frequencies,  signal-arrival  angles,  and  sound 
speed.  Sensitivity  of  the  tracking  algorithm  to  these  types  of  measurement 
errors  is  studied  by  considering  a  variety  of  S-R  scenarios.  In  all  config¬ 
urations,  small  perturbations  of  all  algorithm  inputs  produce  small  variations 
of  output  estimates.  For  example,  in  one  scenario,  percentage  depth,  range, 
and  speed  errors  were  less  than  8%  for  relative  errors  in  periods  of  up  to 
20%.  However,  some  outputs  are  shown  to  be  more  sensitive  to  larger  varia¬ 
tions  in  certain  inputs.  Further  insight  into  algorithm  robustness  is  gained 
by  considering  special  tracking  scenarios,  such  as  motion  confined  to  a 
vertical  plane,  in  which  fewer  than  six  source  parameters  are  unknown.  Our 
sensitivity  analysis  for  these  situations  corroborates  the  algorithm  robust¬ 
ness  for  input  variations  of  reasonable  magnitude. 

For  the  special  tracking  cases,  we  develop  tractable  formulas  for  output 
errors  for  non-random  and  random  input  perturbations.  These  formulas  are  used 
to  obtain  both  qualitative  and  quantitative  information  concerning  relative 
effects  of  various  descriptor  errors.  For  the  random  variations,  we  analyze 
effects  of  phase-roll  period  and  received-frequency  perturbations  on  relative 
errors  in  source  depth  and  speed.  The  influence  of  the  correlation  between 
perturbations  in  these  descriptors,  which  are  both  related  to  travel-time 
variations  due  to  S-R  motion,  on  the  correlation  between  relative  output 
errors  is  determined.  For  example,  if  these  inputs  are  perfectly  correlated, 
then  the  same  result  for  correlation  holds  for  the  outputs.  If  the  standard 
deviation  for  phase-roll  period  variations  is  much  greater  than  that  for 
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frequency  perturbations,  the  degree  of  correlation  between  output  errors  is 
nearly  identical  to  the  input-error  correlation* 

We  also  investigate  the  influence  of  a  specific  corrupting  mechanism, 
namely  Gaussian  white  noise,  on  the  measurement  of  period  and  received 
frequency*  Cramer-Rao  lower  bounds  for  variances  of  unbiased  estimators  for 
relative  errors  in  source  depth  and  speed,  in  terms  of  signal-to-noise  ratio, 
are  developed.  These  bounds  on  tracker  performance  are  derived  by  using  a 
two-stage  estimation  procedure  applied  to  our  algorithm.  Examples  for  special 
S-R  scenarios,  which  typically  involve  very  small  bounds  on  speed  errors  and 
percentage  depth  errors  less  than  15%  for  moderate  to  high  signal-to-noise 
ratios,  indicate  good  tracker  performance. 
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TABLE  I.  Values  for  parameters,  descriptors,  and  outputs  for  scenarios 
(A)-(E)  (c  =  1520  ms-1,  At2  =  120  s). 


Parameter, 
descriptor,  or 
output. 

(A) 

(B) 

Scenario 

(C) 

(D) 

(E) 

Dr  (m) 

200 

200 

150 

200 

200 

vr (ms"1 ) 

5 

5 

10 

10 

10 

4*1  (deg) 

45.0 

100.0 

45.2 

1  .4 

2.9 

^2 

31  .4 

124.2 

106.1 

2.6 

4.1 

(s /cycle ) 

22.0 

6.2 

12.4 

27.4 

21.3 

p2 

.  9.0 

11.4 

22.1 

8.5 

10.9 

fr1 (Hz) 

503.49 

994.40 

503.74 

504.93 

501 .64 

fr2 

503.11 

991 .82 

499.06 

504.93 

501 .64 

Ds  (m) 

300 

300 

300 

300 

300 

R-l  (m) 

3000 

2000 

2000 

4000 

2000 

*2 

1779 

3316 

131 1 

2200 

1400 

vs (ms-1 ) 

10 

10 

5 

5 

5 

3r1 (deg) 

45.0 

100.0 

-45.0 

0 

0 

Pr2 

31 .2 

124.3 

-106.2 

0 

0 

8sl (deg) 

-45.0 

140.0 

-30.0 

0 

180.0 

&s2 

-58.8 

164.3 

-91.2 

0 

180.0 

fg (HZ ) 

500 

1000 

500 

500 

500 
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(a)  Vertical  projection  at  time  t^  of  S  and  R  motions  and  ray 

arrivals.  (b)  Horizontal  projection  of  S  and  R  tracks  between  times 
tj  and  t2«  (c )  Angles  associated  with  direct  ray  at  R. 

Top  view  (time  t^)  of  tracking  scenarios  (A)-(E). 

Percentage  depth,  range,  and  speed  errors  versus  phase-roll  period 
perturbation  (s/cycle)  at  time  t2*  (a)  scenario  (A),  (b) 

scenario  (B) .  Received  frequency  perturbation  =  0.1  Hz. 

For  scenario  (C)  at  time  t2 ,  percentage  depth,  range,  and  speed 
errors  versus  (a)  phase-roll  period  perturbation  (s/cycle)  (received 
frequency  perturbation  0.1  Hz),  (b)  received  angle  perturbation 
(deg) . 

Bearing  and  course  angle  errors  (deg)  at  time  t2  versus  received 
angle  perturbation  (deg),  (a)  scenario  (Bf,  (b)  scenario  (C) . 
Percentage  depth,  range,  and  speed  errors  at  times  t^  and  t2  versus 
phase-roll  period  perturbation  (s/cycle),  (a)  restricted  scenario 
(A) ,  (b)  scenario  (D) .  Received  frequency  perturbation  0.1  Hz. 

Percentage  depth,  range,  and  speed  errors  at  times  and  t2  versus 
received  frequency  perturbation  (Hz),  (a)  restricted  scenario  (A) 0 

(b)  scenario  (E) .  Phase-roll  period  perturbation  1  s/cycle. 
Correlation  coefficient  of  relative  depth/speed  errors  versus  corre~ 
lation  coefficient  of  period/received  frequency  perturbations,  for 
several  standard-deviation  ratios  ic  in  restricted  scenarios. 
Cramer-Rao  lower  bounds  for  percentage  normalized  standard  deviation 
of  depth  and  speed  versus  signal-to-noise  ratio  (dB) . 
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